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Abstract 


At the heart of tj Carinae’s spectacular “Homunculus” nebula lies an extremely luminous (i Total > 5 X 10 6 L Q ) colliding wind binary with a highly 
eccentric ( e ~ 0.9), 5.54-year orbit (Figure 1). The primary of the system, a Luminous Blue Variable (LBV), is our closest ( D ~ 2.3 kpc) and best 
example of a pre-hypernova or pre-gamma ray burst environment. The remarkably consistent and periodic RXTE X-ray light curve surprisingly 
showed a major change during the system’s last periastron in 2009, with the X-ray minimum being ~50% shorter than the minima of the previous 
two cycles 1 . Between 1998 and 2011, the strengths of various broad stellar wind emi s sion lines (e.g. Ha, Fe II) in line-of-sight (l.o.s.) also decreased 
by factors of 1.5 - 3 relative to the continuum 2 . The current interpretation for these changes is that they are due to a gradual factor of 2 - 4 drop in 
the primary’s mass-loss rate over the last ~15 years 1,2 . However, while a secular change is seen for a direct view of the central source, little to no 
change is seen in profiles at high stellar latitudes or reflected off of the dense, circumbinary material known as the “Weigelt blobs” 2, 3 . Moreover, 
model spectra generated with CMFGEN predict that a factor of 2 - 4 drop in the primary’s mass-loss rate should lead to huge changes in the observed 
spectrum, which thus far have not been seen. Here we present results from large- (±1620 AU) and small- (±162 AU) domain, full 3D smoothed 
particle hydrodynamics (SPH) simulations of tj Car’s massive binary colliding winds for three different primary-star mass-loss rates (2.4, 4.8, and 
8.5 X 10" 4 M Q /yr). The goal is to investigate how the mass-loss rate affects the 3D geometry and dynamics of i; Car's optically-thick wind and 
spatially-extended wind-wind collision (WWC) regions, both of which are known sources of observed X-ray, optical, UV, and near-IR emission and 
absorption. We use two domain sizes in order to better understand how the primary’s mass-loss rate influences the various observables that form at 
different length scales. The 3D simulations provide information important for helping constrain t] Car's recent mass-loss history and future state. 
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The 3D SPH Simulations 


We use an improved version of the 3D SPH code in 
[4, 5]. Radiative cooling and radiative forces, 
including inhibition, are now implemented. The 
stellar winds are parameterized using the standard 
beta-velocity law v(r) = v QO (l - RJr)P, with /? = 1, v„ 
the wind terminal speed, and R+ the stellar radius. 
The binary orbit is set in the xy plane, with the 
origin at the system center-of-mass and the orbital 
major axis along the x-axis. The outer simulation 
boundary is set at either r = ±10a or r = ±100a. The 
adopted model parameters (Table 1) are consistent 
with those derived from the observations. By 
convention, <f> = 0, 1, 2, ..., is defined as periastron 
and <p = 0.5, 1.5, 2.5, ..., apastron. 


Figure 2 (right): Top row: Slices from the 3D simulations 
at a phase near apastron showing log density in the orbital 
xv plane for 3 assumed primary mass-loss rates. The stars, 
WWC cavity, and Weigelt blobs are marked. The dashed 
line and eye show the observer’s projected l.o.s. Axis tick 
marks correspond to an increment of 10a = 1 54 AU ~ 0.07”. 
Bottom row: Same as top row, but in the xz plane. 

Far right panel: Illustration of t] Car’s orbit (upper left 
inset, yellow) on the sky relative to the Homunculus for the 
binary orientation derived by Madura et al. (2012) (i = 138°, 
co = 263°, and PA Z = 317°). The red, green, and blue arrows 
indicate the orbital major (+x), minor (+y), and angular- 
momentum (+z) axes, respectively. North is up. 


Orbit projected on sky relative to Homunculus 


The Spatially-Extended, Time-Variable 
Interacting Wind Structures 


Collapse of the Inner Wind-Wind Collision Zone 
During Periastron Passage 


-a.S k 1G" M. ftp 




Figure 4: Slices showing density, temperature, and wind speed on a logarithmic scale (columns, left to right, cgs units) in the orbital xy plane at 
seven phases around periastron (rows) taken from the r = ±10a SPH simulations of rj Car assuming primary mass-loss rates of 8.5 x 10 -4 M Q /yr 
(left three columns) and 2.4 x 10 -4 M Q /yr (right three columns). All plots show the inner ±la region. Axis tick marks correspond to an increment 
of 0. 1 a ~ 1 .55 AU. The orbital motion of the stars is counterclockwise. The LBV primary is to the left and the companion to the right at apastron. 
In each simulation there is a “collapse” of the WWC zone between the stars. The hottest gas near the WWC apex vanishes during this time. The 
higher the primary mass-loss rate, the sooner the collapse occurs before periastron, and the later the recovery of the hot gas after periastron. The 
collapse appears to be due to the increased effectiveness of radiative cooling in the post-shock wind of the companion at phases around periastron. 


Figure 3: Slices showing log density (left two columns) and log temperature (right two columns) at apastron (0 
periastron (<f> = 3.0) in the xy (top two rows), xz (middle two rows), and yz (bottom two rows) planes taken from t 
3D SPH simulations of t] Car assuming primary mass-loss rates of 8.5 x 10 -4 M Q /yr and 2.4 x 10~ 4 M Q /yr (colun 
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ary Compani 


Mass (M 0 ) 

90 

30 

Radius (R 0 ) 

60 

30 

Mass-loss Rate (10 -4 M G /yr) 

8.5, 4.8, 8t 2.4 

0.14 

Wind Terminal Speed (km/s) 

420 

3000 

Wind Momentum Ratio /] 

0.12, 0.21, & 0.42 

Orbital Eccentricity e 

0.9 


Orbital Period 

2024 days 

Semimajor axis length a 

15.45 AU 
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